General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



3 (H1SI-TH-85258) ULTBAflOLIT SEFCTFCSCCFY OF H03- 20875 

THI EH1GHTEST SOPEBGI AHTS IS B 3 1 AFE H33 
Pinal Technical Report (Minnesota Osiv.) 

«HI p HC A03/HF A0 1 CSCL 03A Unclas 

G3/89 03002 


FINAL TECHNICAL REPORT 
NASA/NAG 5-191 


Ultraviolet Spectroscopy of the Brightest Supergiants in M31 and M33 

by 

R. M. Humphreys, T. R. Gull, and S. D'Odorico 


University of Minnesota 
Department of Astronomy 
116 Church St. S. E. 
Minneapolis, Minnesota 55455 






IUE and Groundbased Observations of the Hubble-Sandage 
Variables In M31 and M33 


R.M. Humphreys^ and C. Blaha 
University of Minnesota 

S. D’Odorlco 

European Southern Observatory 

T.R. Gull 

NASA-Goddard Space Flight Center 
and 

P. Benvenuti 

Astronomy Division, European Space Agency 


1 


Guest Observer, Mt. Palomar Observatory. 


Abstract 


Ultraviolet spectroscopy from the IUE, In combination with groundbased 
visual and infrared photometry, allows us to determine the energy distributions of 
the luminous blue variables, the Hubble-Sandage variables, in M31 and M3 3 . The 
observed energy distributions, especially In the ultraviolet, show that these 
stars are suffering interstellar reddening. When corrected for interstellar 
extinction, the integrated energy distributions yield the total luminosities and 
black body temperatures of the stars. The resulting bolometric magnitudes and 
temperatures confirm that these peculiar stars are indeed very luminous, hot 
stars. They occupy the same regions of the Mbol vs. log T e diagram as do n Car, 

P Cyg and S Dor in our galaxy and the LMC. 

Many of the Hubble-Sandage variables have excess infrared radiation which is 
attributed to free-free emission from their extended atmospheres. Rough mass 
loss estimates from the infrared excess yield rates of 10~3 M^/yr. The ultraviolet 
spectra of the H-S variables are also compared with similar spectra of n Car, P Cyg 
and S Dor. 


Subject headings: galaxies: individual - stars: supergiants - stars: luminosities - 


stars: variables. 
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lit INTRODUCT ION 

The Hubble-Sandage (H-S) variables are among the most luminous stars 
known. These blue variables are spectroscopically and photometrically similar 
to the highly luminous, peculiar stars n Car and P Cyg in our galaxy and S Dor 
in the Large Magellanic Cloud. They are included in the class of stars termed 
"S Doradus variables" by Kukarkin et_al. (1974). Five blue variables in M31 
and M33 were the original H-S variables (Hubble and Sandage 1953): 

AF And (Var 19) in M31 and variables A, B t C and 2 in M33. Other luminous, 
blue stars which resemble the original group are also included as H-S variables; 
Hubble Var 15 and Var A-l (Rosino and Bianchini 1973) and AE And (Sandage and 
Tamaann 1974) in M31 and Var 83 in M33 (van den Bergh, Herbst and Kowal 1975). 

These variable stars share common spectroscopic and photometric charac- 
teristics. They display irregular variability consisting of extended maximum 
and minimum phases, frequently lasting several years, and separated by rela- 
tively short transition phases. Their light curves may be found in Hubble and 
Sandage (1953) and Rosino and Bianchini (1973). The spectral characteristics 
of the H-S variables are discussed by Humphreys (1975, 1978). These variables 
display a strong, hot continuum extending into the ultraviolet, no apparent 
Balmer discontinuity and strong emission lines of H, He I and Fe II 

and [Fe II]. They show ultraviolet excess radiation and some have Infrared 
excesses (Humphreys and Warner 1978). In AF And, Var A-l, Var 83, and Var 2, 
the infrared excess may be due to free-free emission from extended 
atmospheres, while in Var A it is very likely caused by circumstellar dust. 

Current models suggest that the H-S variables are hot, massive supergiants 
surrounded by circumstellar envelopes (Humphreys 1978, Gallagher, Kenyon and 
Hege 1981, Wolf, Appenzeller and Cassatella 1980). Their photometric behavior 
combined with their low surface gravities, suggest that their atmospheres may 
be unstable in a manner similar to j) Car (Davidson 1971). A better understanding 



-3- 


•ORIOHAL PAGE » 
Of POOR QUALITY 


of their evolutionary status and their role in massive star evolution requires 
information on their total luminosities and temperatures. 

For this reason ultraviolet spectra were obtained with the International 
Ultraviolet Explorer (IUE) of the brightest H-S variables. The ultraviolet fluxes in 
combination with visual and infrared photometry* allow us to determine their 
energy distributions, luminosities and temperatures. The ultraviolet spectra 
can also be compared with possibly related objects like n Car, S Dor, and F Cyg. 

The observations and the resulting energy distributions are described in sec- 
tions 2 and 3. The ultraviolet spectra are presented in section 4. In section 5 
the importance of the Hubble-Sandage variables in massive star evolution is 
discussed. 


II. OBSERVATIONS 

Although the H-S variables are blue and have ultraviolet excesses, they 
are typically 16 to 17 magnitude in the visual, and only the very brightest 
can be observed with the IUE. Only five of the H-S variables were bright 
enough and then only marginally. These stars were all brighter than V - 17 mag 
and most were near V » 16 mag at the time of observation. Long (LWR - 1800-3255 A 
and short (SWP «* 1135-2085 A) wavelength low dispersion spectra were obtained of 
AF And in M31 and Var 83 and Var 2 in M33, but only long wavelength spectra were 
observed for Var B in M3 3 and Var A-l in M31. Because of their faintness these 
stars could not be seen with the acquisition telescope and had to be acquired 
by blind offset from nearby stars which were typically V s 12 to 14 mag. Table 1 
is a journal of the ultraviolet observations, and includes the 1950 coordinates 
for these objects. The positions were measured from plates of M31 and M33 at 
Kitt Peak and ESO, and are based on coordinates of SAO stars. The accuracy of 
these coordinates is 1 arc sec from both the Internal accuracy of 

the measurements and a comparison of the two independent determinations. 


f 
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All of the 1981 IUE spectra ware obtained at NASA's Goddard Space Flight 
Center, and were reduced according to standard procedures. Known reseau Marks 
and hot pixels were removed from the spectra and patched over before any Integrated 
fluxes were measured. The 1982 IUE spectra of Var 83 and AF And were recorded 
at the Vlllafranco Satellite Tracking Station of the European Space Agency. 

The geometrically and photometrically two-diaanslonal corrected Images were 
converted to the format of the ESO two-dimensional data reduction system (IKAP) 
and analysed In Garching. 

Broadband fluxes were measured from the calibrated ultraviolet spectra 
for three passbanda - a short wavelength (SW), a long wavelength (IM), and one 
(MW) centered on the 2200X feature. The central wavelength and width of each 
passband are given below 



X (y) 

AX(y) 

SW 

0.160 

0.1225-0.1975 

MW 

0.220 

0.2050-0.2350 

LW 

0.270 

0.2350-0.3050 


and the corresponding fluxes for the five variables are in Table 2. Narrower 
passbands were not used because of the stars' faintness. 

All cf the available visual and infrared photometry is summarised in Table 3. 
The 1976 and 1977 observations are from Humphreys (1978) and Humphreys and Warner 
(1978). Additional UBVR photometry was measured in I960, 1981 and Aug. 1982 
with the Mark I 'computer photometer' on the KPNO 2.1-meter talescope and in 
Oct. 1982 with the pulse-counting photometer on the UM-UCSD 1.5-meter on Mt. 
Lemmon. Both visual photometers are aperture chopping Instruments with simul- 
taneous sky measurements. The . JHK observations were made with the 

lndium-antimonide system on the Mt. Pslomar 5-meter in Oct. 1980. 

III. THE SPECTRAL ENERGY DISTRIBUTIONS 

Fluxes ranging from the ultraviolet to 2.2y are now available for five of the 
H-S variables, and It is possible to discuss their energy distributions snd a 
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determination of their bolometrlc luminosities and temperatures. However, the 

observations must first be corrected for interstellar extinction, both 

foreground and Internal reddening. These stars are in young stellar associations 

in the spiral arms of M31 and M3 3, and It is very likely that they suffer some 

reddening due to dust Inside their galaxies. Indeed their observed energy 

distributions show that these stars are suffering some reddening. 

Because of their peculiar emission line spectra, we cannot use the 

observed multicolor photometry to estimate their reddening. An alternative 

method is available from the neutral hydrogen maps which provide a means of 

estimating the visual extinction along the line of sight. Knapp, Kerr and Hose 

(1973) and Savage and Jenkins (1972) have shown that the column density of 

21 

neutral hydrogen in our galaxy is proportional to ■ 5.6 x 10 Eg_ v ). 

The n _ is then determined along the line of sight to each star from the HI 
HI 

maps published by Roberts (1966) for M31 and by Wright, Warner and Baldwin (1972) 
for M33. The corresponding visual extinction, A^» is found from the above 
relationship assuming that the ratio of total to selective extinction is 3. 

We do not know the exact location of the stars with respect to the neutral hydrogen 
along the line of sight and it is rather unlikely that all of the extinction 
estimated from the HI column density should be applied; consequently, we have 
assumed that on the average the stars experience only half the extinction 
from the neutral hydrogen. We also know thjt the minimum A^ due to foreground 
reddening is 0.4 mag in front of M31 (van den Bergh 1968) and 0.3 mag in front 
of M33 (van den Bergh 1968, Humphreys 1980). The adopted A v is then defined 
as the foreground A and half that from the neutral hydrogen column density. 

This information is included in Table 4, together with the visual magnitude 
closest to the time of the ultraviolet observations and the resulting absolute 
visual luminosity from the adopted true distance moduli of 24.1 mag (van den Bergh 
1976) and 23.9 mag (Humphreys 1930) for M31 and M33, respectively. 

The broadband energy distributions from the ultraviolet to 2.2y are 
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shown In Figures 1 through 6 for the five H-S variables with IUE observation*. 

The visual and infrared photometry closest to the time of the ultraviolet 
data are used for the energy distributions. This presents a small problem 
for some of the stars which were varying, because it was not always possible 
to have nearly simultaneous observations, but the colors of most of the stars 
show only small fluctuations even though the visual magnitude may be varying 
somewhat. Only one complete set of Infrared photometry (JHK) is presently 
available for moat of the stars, and these same Infrared colors were assumed 
to be applicable to visual photometry obtained later. The original observations 
and the data corrected for the visual extinction defined above with the galactic 
ultraviolet extinction curve from Nandy et al . (1975) are both shown in the plots 
of the energy distributions. A blackbody is then fit through 

the dereddened observations using the galactic extinction curve. With these 
adopted blackbody curves the energy distributions were Integrated from 0.09u 
to 2. 2y to determine the total flux of the star and consequently its total 
luminosity. The blackbody temperatures and the luminosities (M^^) determined 
in this way are also included in Table 4. 

We have investigated the effects of varying the adopted extinction on the 
blackbody temperatures and luminosities. It is obvious from the observations 
that all of the stars require some correction for interstellar extinction. The 
minimum extinction is that due to foreground reddening. A blackbody curve fit 
to these stars' observations when corrected for minimum extinction can usually 
be matched reasonably well with the ultraviolet photometry but gives a poor 
fit to the red and infrared observations, often with a very large infrared 
excess. (Evidence for free-free emission from some of the 

H-S variables is discussed below.) Figures 3 and 8 for Var A-l 

and Var 15, respectively, show the effects of varying the extinction between 
the adopted value defined above and the maximum value from the neutral hydrogen 
colimin density. For Var A-l higher extinction is probably appropriate 
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and results In a batter blackbody fit to the observations. However, for nost 
of the H-S variables the maximum extinction gives a much poo*, or blackbody fit 
to the observations, often resulting in a very large depression of the visual, 
red and infrared photometry below the blackbody curve. Thus it seems that our 
adopted extinction based on the foreground reddening and the average extinction 
from the neutral hydrogen column density is justified for most of the stars. 

Figure 1 for AF And illustrates the effects of varying the ultraviolet 
extinction curve. It is clear that the galactic extinction curve over corrects 
the flux in the 2200A feature. This is the case for all of the H-S variables in 
both M31 and M3 3 and the excess may even be a little higher for stars in M33. On the 
average the 2 2 00 A feature from hot stars in the LMC is somewhat lower, and 
when the ultraviolet extinction cu*.» -i for the LMC (Handy et al . 1980) is used 
the "excess" at 2200l is largely eliminated for most of the stars, although it 
is still slightly above the blackbody curve for stars in M33. Use of the IMC 
extinction curve also results in higher temperature blackbody fits to the 
observations and consequently higher luminosities. It is perhaps premature 
to draw conclusions about the ultraviolet interstellar extinction in M31 and 
M33 on the basis of these data because the 2200& region has low sensitivity 
vith the IUE detector and the calibration at these low flux levels is uncertain. 

This is certainly a problem for future work when higher quality data are available 
for stars in more distant galaxies. 

When the best-fit blackbody curve is placed through the ultraviolet and 
visual observations for some of the H-S variables, it is readily apparent that 
their infrared fluxes lie somewhat above the blackbody curves. It is very 
likely that these stars have some excess infrared radiation due to free-free 
emission from their extended atmospheres. The low resolution IUE spectra are 
Inadequate to comment on evidence for mass loss from the line profiles. However, 
a rough mass loss (M) estimate can be determined from the J-K color excess due 
to free-free emission from Hyland's (1979) empirical calibration of E^(J-K) 
versus M/V^ where V m is the terminal velocity. The adopted intrinsic J-K 
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color* are from Johnson (196' using the best-fit blackbody temperature from 
the energy distributions, to approximate the true stellar temperature. Because 
we cannot measure directly for these stars we adopted the value for P Cyg 
of 300 km/sec from Barlow and Cohen (1977). This assumption is reasonable 
because the H-S variables are spectroscopically and photometrically similar 
to P Cyg (see section V). R81 In the LMC, another member of this class of 

variables, has a V m of 250 km/sec (Wolf et al . 1981). The resulting mass 
loss rates are given in Table 4 and for all the stars are around 10 ^ M^/yr 
which Is comparable to the mass loss rates for the ocher luminous variables 
P Cyg, R81, R71 and S Dor. Although these rates are quite reasonable, they 
must be considered preliminary until higher resolution ultraviolet spectra 
become available. 

The energy distributions for the Individual stars are discussed below. 

M31 - AF And - An LWR spectrum was obtained with IUE In November 1981, and 

In Oct. /Nov. 1982 both SWP and LWR observations ware successfully made. The 

latter observations are used for the energy distribution in Figure 2. it is 

fortunate that there is nearly simultaneous visual photometry from Aug. and 

Oct. 1982, although it is necessary to use the Infrared colors fro*.’ the 1980 

observations. A 2^ 000 K blackbody gives a good fit to the ultraviolet and 

visual data; however, the Infrared flux is significantly above any blackbody 

fit to the short wavelength observations. The same la true with the 1980 and 

1981 observations. This Infrared excess radiation, attributed to free-free 

emission in the star's extended atmosphere corresponds to a mass loss rate of 
-5 


-9- 


OVHQBIAL PAQK VS 
OF POCK QUALITY 


Var A-l - Only ’long wavelength* (LVR) observations from IUE arc 
available for Var A-l and its energy distribution is shown in Figures 3s end 3b. 
The visual photometry from Dec. 1981 was obtained only one month after the 
ultraviolet spectrua, As with AP And it is difficult to fit the infrared 
observations with the visual and ultvavlolet data. Figure 3s shows a 15,800 K 
blackbody fit to the data corrected for the adopted extinction. However 
Gallagher, Kenyon and Hege (1981) suggest that Var. A-l say be In a region of high 
extinction. For this reason, the observations are also shown in Figure 3b 
corrected for extinction based on the aaxiaua value fro* the neutral hydrogen 
coluan density. A 27,300 K blackbody then gives a tighter fit to the 
observations at all wavelengths with a saall excess in the Infrared . A axes 
loss rate of 1.2 x 10 ^ M^/yr ia derived from this free-free emission in the 
infrared. 

M33 - Var 83 - This scar seems to be the hottest in the group at this time. 
It is also now the visually brightest; consequently, it has the largest nuaber 
of IUE observations. It is worth noting that Var 83 has been slowly getting 
brighter over the past six years and it is now the brightest it has ever been 
observed (V - 15.4 aag). Both LWR and SUP observations of Var 83 were obtained 
in June and November 1981 and in Oct. /Nov. 1982. During this tlae the star 
brightened in the visual by nearly a magnitude. This aakes it possible to 
compare the star's temperature and total luminosity as it brightens. 
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Th« energy distribution* era shown in Figures 4a end 4b for the Nov. 1981 
end Nov. 1982 observations, respectively. The 1981 vlpual photometry wee 
obtained only on* month later, but the infrared date ere from the previous 
year. The best blackbody fit to the 1981 observations gives a temperature 
of 29100 K with an infrared excess, most likely the result of free-free 
amission. The corresponding mass loss rate le 3 x 10~~ M # /yr. The 1982 
observations show a significant change in the fluxes for both ultraviolet 
passbands. The ehort wavelength (.16w) flux decreased while whs star 
brightened at the longer wavelength (.27y). The energ$ distribution clearly 
shows a shift in the peak flux to longer wavelengths, and a 20,000 K blackbody 
gives a good fit to the 1982 data. As the star brightened in the visual, its 
temperature decreased significantly by nearly 50Z, however the total luminosity 
of Vsr 83 remained essentially the asms; the bolometric magnitudes are 
•11,3 mag and -11,1 mag in 1981 and 1982, respectively. 

Var 2 - Both LWR end SWF spectra were obtained with IUE in June 1981, but 
Var 2 decreased in brightness by nearly one magnitude between Nov. 1980 and 
Dec. 1981. Figure 3 shows the energy distribution with both sets of visual 
photometry. The ultraviolet data, taken in June 1981, fit quite well with 
the 1980 visual and infrared data with a blackbody temperature of 22,800*K. 

The Dec. 1981 photometry appears too low to fit the ultraviolet dace. Ther- ore 
it is suspected that moat of the decline in brightness occurred between June 
and Dec. 1981. Indeed, if Var 2 had been as dim in June as the December 
observations indicate, it could not have been observed with IUE. 

The infrared excess yields . mass loss rat* of about 3.5 x 10” 5 Mg/yr. 

Var B - Only the LWR spectrum is available from IUE, and no infrared 
observations were made of Var B in 1980. With this limited wavelength 
coverage only a rough temperature estimate can be obtained; however, the 
available photometry definitely Indicates a hot blackbc-dy. A 21,900*K 
blackbody give a good fit to the data in Figure 6. 
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The remaining four H-S variables were Coo faint at the tine for observation 
with IUE. However for thu purpose of completeness It la possible to make 
sone estlnate of the temperature and luminosity of AE And and Var 15 in M31 
from the visual and Infrared photometry. Var C In M33 has been omitted from 
the discussion in this paper because of its limited wavelength coverage; only 
visual photometry is available. Var A is the most unique and potentially, the 
most interesting of the H-S variables. It is the only one which apparently has 
a clrcumstellar dust shell. Its peculiar energy distribution is discussed here, 
but because of its very high reddening it is not possible to determine the ' 
temperature and luminosity without further long wavelength observations. 

M31 - AE And - This blue variable has photometric data from the visual 
through the infrared. The 1980, 1981 and 1982 visual observations are all 
shown in Figure 7 with the 1980 infrared photometry. The observed energy 
distribution is fit reasc-*bly well by a 15,000 K blackbody. It is possible 
that a higher temperature blackbody could be drawn through the visible 
points, with tli« infrared partially attributed to free-free Mission. Owiousiy, 
ultraviolet data la needed for a more accurate determination of the 
temperature of AE And. 

Var 15 - This star appears to be the coolest of the H-S variables. The 
energy distribution in Figure 8* is best fit by a blackbody of 8300 K. Evan 
with the maximum visual extinction from the HI data (Figure 8b) applied, a 
relatively low temperature is still indicated, 12,400 K. Since a low temperature 
is unusual for this group of iters, the 1973 and 1977 observations were also 
examined and they also ‘suggest a temperature near 8000 K. 
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H33 - Var A - This peculiar star is the most Interesting member of the 

H-S group. It may veiy well be an extragalactic n Car. Its light curve 

(Hubble and Sandage. 1953, Rosino and Bianchini 1973) Is quite similar to n Car. 

In late 1950, Var A reached maximum brightness (M * 15.7 mag) with a color 

PIS 

index of +0.33 mag. After a very rapid decline followed by a secondary peak, 

the color index had reddened to +1.5 and the star had dimmed to M * 19 mag. 

pg 

It has remained faint ever since (see Table 3 ). Humphreys and Warner (1978) 
found that Var A has a large infrared excess with a V-K of 4.8 mag. Var A's 
photometric behavior plus its Infrared radiation suggests the ejection of one 
or more shells in 1950-53, followed by the formation of drcumsteilar dust. 

Figure 9 shows its energy distribution. Clearly the infrared 
flux is increasing with wavelength and is highly suggestive of circumstellar 
dust. The observations corrected for Interstellar extinction determined as 
for the other variables are shown, and it is obvious that Var A may have 
additional reddening perhaps due to circumstellar dust. To test this idea 
a value for circumstellar extinction is estimated. Assuming dust formed 
after Var A reached maximum, much of the presently observed red color may 
be attributed to circumstellar dust. The color index at maximum, transformed 
to B-V, is assumed to be due to the star plus some interstellar reddening, 
and the extinction from circumstellar dust is defined as: 

(I Wcs ■ ( »- V) nov - (B - V) max ’ 

A B-V of 1.01 mag is adopted from the observations in Table 3 and the 

excess attributed to circumstellar dust is .57 mag. Assuming the same ratio 

of total ; o selective extinction (R) as for interstellar dust yields A of 

CS 

1.71 mag. Th« value of R may be larger than 3 for the circum- 

stellar dust, but it is at present unknown. This then gives a total reddening 

of 2.73 mag (A + A ). The energy distribution corrected for this total 
V IS V CS 
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reddening i also shown in Figure 9. These results ere of course ell very 
preliminary, but it is very interesting that with a correction for circuastellar 
reddening Var A is still brightening at K(2.2y). Certainly longer wavelength 
observations especially at 10y are needed. 

IV. THE ULTRAVIOLET SPECTRA OF THE HUBBLE-SANDAGE VARIABLES 

Several of the brightest members of the 'S Doradus* class of variables 
have been observed previously with 1UE. Wolf et al. (1980, 1981a, b) have discussed 
the ultraviolet spectra of S Dor, R71, and R81 in the LMC.. Cassatella et al. (1979) 
Luud and Sapar (1980), and Laniers et al. (1983) have described the ultraviolet 
spectrum of P Cyg, and Cassatella, Giangrande and Vlotti (1979) and Viotti et al. 
(1981) have briefly commented on the spectrum of 0 Car. 

All of these stars are in the Galaxy or the LMC, and are therefore sufficiently 
bright to be observed with the IUE high dispersion camera. Most of the previous 
work, mentioned above, is primarily concerned with an analysis of the line 
profiles and a determination of the mass loss rates. There has been no previous 
attempt to compare the energy distributions and the main spectral features in the 
ultraviolet region as derived from the low dispersion data. Therefore, we have 
also obtained low dispersion spectra of u Car, P Cyg, and S Dor for a comparative 
discussion with our ultraviolet spectra of the H-S variables . for which only a few 
features can be recognized due to the low signal to noise. The combined IUE 
short and long wavelength spectra of n Car, P Cyg, and S Dor are shown in 
Figure 10a with the principal line Identifications from the above references. 

The ultraviolet spectra of AF And in M31 and Var83 in M33 from two different 
times are shown in Figure 10b. 

The observed energy distributions of P Cyg and 0 Car are strongly affected 
by interstellar reddening and both stars showed a marked 2200 A depression. Both 
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th« short and long vavalsngth spactra of n Car show strong features, many In 

absorption, which ara bland® of aatallic linaa and ara difficult to dlsantangla 

at low dispersion. Tha spectrum of r| Car shows a strong MgXI doublet at 2793 A 

and 2803 A, a vary broad absorption cantarad at 2600 X due mainly to Veil 

multlplat 1, and a sharp drop balow 2415 A, attributed to a combination of 

absorption by Fell multlplat 2 and extinction by tha 2200 A feature.. In tha short 

wavelength spectrum, CIV is not as prominent as three strong absorption features. 

• • • 

centered at 1565 A, 1620 A, and 1700 A to which tha main contributors ara probably 
low ionization ions of Fa, Mg, Mi and Al. 

In P Cyg, tha absorption lines in tha 2000-3000 X region ara not as prominent 
as in n Car, only the Mgll doublet and tha Fell multlplat 1 are readily apparent. 
The same is true for S Dor where only Mgll is strong, and tha fainter lines have 
been identified with Fall by Wolf at al. (1980). The same authors show that in 
the short wavelength range CIV at 1500 A is the strongest feature. 

The change in the slope of the continuum of Var83 between June 1981 and 
November 1982, discussed in tha previous section, is easily recognised in the 
ultraviolet spectra shown in Figure 10b. However, only a few spectral faatures 
can be definitely identified. At short wavelengths. Sill + SII at 1260 A, 

Si II + 01 at 1303 A and CII at 1335 A are probably present and are mainly 
interstellar in origin. SUV 1394-1403 A and CIV 1550 A are not clearly 
Identified. At long wavelengths three absorption features are seen centered 
at approximately 2400 A (Fell, 2), 2600A (Fell ,1) and 2800 A (Mgll). In this 
respect, the spectrum of Var 83 is more like that of n Car than S Dor. In 
AF And we cannot clearly Identify any features. 
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V. DISCUSSION - THE HUBBLE-SAHDAGE VARIABLES AND MASSIVE STAR EVOLUTION 

The IUE date presented in this paper conflm the high luminosities and 
high temperatures suspected from the earlier visual spectra and photometry of 
the H-S variables. Their bolometric luminosities and absolute visual luminosities 
corresponding to the time of the ultraviolet observations are summarised in Table 
4 together with the blackbody temperatures from their energy distributions. Hy 
ranges from -7.0 to -9.6 mag and from -8 to ■ -11.5 mag. The brightest 
magnitudes recorded in the literature (see Humphreys 1978 for references and 
the data) are used to derive the visual luminosities at 'maximum light' which 
range from -9.0 mag to perhaps as bright as -11 mag. Of course some of these 
stars may not have heen observed at their maximum brightness and it is 
unfortunate that we do not know their corresponding bolometric luminosities. 

The positions of the H-S variables on the Mg o ^ vs. temperature H-R 
diagram are shown in Figure 11 using the bolometric luminosities and 
blackbody temperatures in Table 4. This H-R diagram from Humphreys (1983) 
is a composite for the brightest stars in six local Group galaxies. Only 
stars brighter than Mgol " -9.0 mag are shown and the evolutionary tracks 
for massive stars are from Maeder (1980). It is clear from their positions 
on the H-R diagram that the H-S variables are among the most luminous stars 
known. Var 83 and Var A~1 lie on or near the empirically determined upper 
luminosity boundary for normal supergiants. They are both in the same 
general region of the diagram where n Car is also found, while AF And, 

Var 2 and Var B have luminosities and temperatures similar to P Cyg. 

Only AE And and Var 15 have relatively low luminosities compared to the other 
H-S variables. Both are presently very faint, perhaps near minimum light, and neither 
has ultraviolet observations. When ultraviolet data are available higher 
temperatures might be indicated. Actually these two stars are rather similar 
to S Dor. and R71 in the LMC (Appenzeller and Wolf 1981) which have luminosities 
between -8.5 and -9.5 mag and temperatures of 9000 to 14000 K. 
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The H-S variables are obviously very massive stars, presumably evolved 

past the core hydrogen burning stage based on their positions on the H-R 

diagram. Because of their irregular variability and probable high mass 

loss rates, we may surmise that they have reached the point in their 

evolution when the expected pulsational instabilities in the core (Schwarzschild 

and Ham 1959, Appenzeller 1970) and high radiation oressure at the surface 

are responsible for rapid mass loss. In the case of n Car, and presumably 

Var A, the mass loss is very likely sporadic and very rapid (Humphreys and 

Davidson 1979) and has resulted in the formation of an extensive dust shell. 

Hone of the H-S variables has quite the high luminosity of n Car (Mg^ a -12 mag), 

but o^ course most of than have only been observed for about 60 years and the 

information on maximum brightness is necessarily Incomplete. Nevertheless, 

P Cyg, S Dor, R71, and the H-S variables appear to be less extreme 

examples of the n Car phenomenon. It is possible that all stars above 

some initial mass, say >60-80 M , pass through an n Car/P Cyg stage during 

o 

which they shed a large amount of matter. 
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Journal of Observations with IUE 


Star 


H 


Image Ho. 

Date 

Exposure Time (min.) 

M31 AF And 

0 h 40 48.6 

+40 

55 

45.0 

LWR 11850 

November 1981 

420 






SWP 18689 

November 1982 

391 






LWR 14756 

November 1982 

316 

M31 Var A-l 

0 h 42 05.6 

+41 

14 

13.8 

LWR 11853 

November 1981 

430 

M33 Var 83 

l h 31 21.7 

+30 19 

16.7 

LWR 10956 

June 193:. 

210 






LWR 10971 

June 1961 

390 






SWP 14342 

June 1981 

165 






LWR 11881 

November 1981 

340 






SWP 15359 

November 1981 

395 






SWP 8247 

October 1982 

380 






LWP 1718 

November 1982 

365 

M3 3 Var 2 

l h 3i 29.0 

+30 

23 

16.5 

LWR 10966 

June 1981 

240 






SWP 14355 

June 1981 

180 

M33 Var B 

l h 31 00.0 

+30 

22 

48.1 

LWR 11861 

November 1981 

360 




Star 

Data 

Flux 
SW(. 16) 

(10 Watte/cm^) 

MW(.22u) LW(.27 ) 

1' 

MIX AF And 

November 1981 

r 

0.80 

2.35 

i 


November 1982 

.2 .40 

0.60 

2.03 


M31 Var A-l 

November 1981 

- 

0.34 

1.14 

* 

M33 Var 83 

June 1981 

5.10 

1.41 

1.50 



November 1981 

4. 84 

1.24 

3.51 



November 1982 

2.93 

1.29 '* 

3.99 


M33 Var 2 

June 1981 

3.74 

1.09 

2.22 


M33 Var B 

November 1981 

— 

1.46 

3.31 



Tabic 3 


OJWWALPAQEIS 
°f POOR QUALITY 


Summary of Visual and Infrared Photometry 


1 


Star 

Date 

V U-B 

(mag) (mag) 

B-V 

(mag) 

V-R 

(mag) 

V-I 

(»ag) 

V-J 

0»ag) 


mm I 

H31 

AE And 

Sept. 1976 

17.00 -0.81 

+0.10 

+0.34 

+0.62 

• • • 

• • • 

a a a 



Oct. 1977 

17.33 -0.78 

+0.06 

+0.24 

* • • 

• a a 

a • • 

a a a 



Nov. 1980 

17.94 ... 

+0.08 

+0.41 

e e • 

e a a 

+0.91 

+0.82 



Dec. 1981 

17.79 -0.85 

0.00 

+0.47 

a a a 

• • • 

a • • 

a a a 



Aug. 1982 

17.57 -0.95 

+0.08 

+0.41 







Oct. 1982 

17.39 -0.96 

+0.05 

a • • 





M31 

Var 15 

Sept. 1976 

17.47 -0.49 

+0.44 

+0.66 

+0.84 

e • a 

a a a 

a a a 



Oct. 1977 

17.69 -0.57 

+0.36 

+0.63 

• • a 

• • a 

a a a 

a ^ a 



Nov. 1980 

17.32 ... 

• • • 

+0.74 

• • • 

+1.41 

+1.59 

+1.96 



Dec. 1981 

17.27 -0.49 

+0.43 

+0.80 

• • • 

• • • 

a a a 

a a a 



Aug. 1982 

17.01 -0.35 

+0.48 

+0.71 







Oct. 1982 

17.35 -0.27 

+0.33 

+0.95 





M31 

AF And 

Sept. 1976 

16.05 -0.85 

+0.11 

+0.43 

+0.74 

• • • 

a a a 

a a a 



Oct. 1977 

16.11 -0.82 

+0.18 

+0.35 

• • • 

• • • 

a a a 

+1.3 



Nov. 1980 

16.42 ... 

-0.01 

+0.26 

a a a 

+0.74 

+1.07 

+1.33 



Dec. 1981 

16.52 -0.92 

+0.13 

+0.37 

• a • 

# • • 

a a a 

a a a 



Aug. 1982 

16.40 -0.90 

40.18 

+0.41 







Oct. 1982 

16.35 -0.91 

40.20 

+0.49 





M31 

Var A-l 

Sept. 1976 

16.26 -0.54 

40.41 

+0.54 

+0.71 

• • * 

a a a 

a a a 



Oct. 1977 

16.25 -0.49 

+0.40 

+0.58 

• • • 

• • • 

a a a 

+1.8 



Nov. 1980 

16.59 

+0.39 

+0.50 

• • • 

+1.13 

+1.36 

+1.47 



Dec. 1981 

16.75 -0.65 

+0.42 

+0.50 

e • • 

• • • 

a a a 

a a a 



Aug. 1982 

16.39 -0.71 

+0.40 

+0.51 







Oct. 1982 

16.60 -1.12: 

40.35 

+0.95: 





M33 

Var B 

Sept. 1976 

17.11 -0.82 

+0.01 

40.93 

• • • 

• • • 

a a a 

a a a 



Oct. 1977 

17.31 -1.02 

-0.11 

+0.05 

• • • 

• • • 

a a a 

a a a 



Dec. 1981 

16.14 -0.75 

+0.22 

+0.37 

• • • 

• • • 

a a a 

a a a 



Aug. 1982 

16.30 -0.39 

+0.14 

+0.34 







Oct. 1982 

15.63 -0.51 

+0.31 

40.53 





M33 Var 2 

Sept. 1976 

18.17 -1.00 

-0.17 

40.46 

• • • 

• • • 

a a a 

a a a 



Oct. 1977 

18.10 -1.03 

-0.19 

+0.42 

• • • 

• • a 

a a a 

a a a 



Nov. 1980 

16.36 ... 

+0.22 

+0.41 

• • • 

+0.59 

+0.89 

+1.05 



Dec. 1981 

17.29 -0.93 

+0.05 

+0.43 

• • • 

a a a 

a a a 

a a a 



Aug. 1982 

17.43 -0.85 

+0.15 

+0.45 







Oct. 1982 

17.63 -1.00: 

+0.01; 

+0.51: 





H33 

Var 83 

Sept. 1976 

16.72 -0.93 

+0.05 

+0.22 

• • • 

• • • 

a a a 

a a a 



Oct. 1977 

16.63 -0.87 

0.00 

+0.30 

• a • 

• • a 

a a a 

+2.5 



Nov. 1980 

16.51 ... 

+0.10 

+0.34 

• • • 

+0.27 

+0.45 

+0.68 



Dec. 1981 

16.16 -0.90 

+0.05 

+0.37 

• • • 

• a • 

a a a 

a a a 



Aug. 1982 

15.67 -0.90 

+0.18 

+0.35 







Oct. 1982 

15.43 -0.83 

+0.20 

+0.24 






&mrnt£mm .... i.i-UffciTi- i ■■ - • - 
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Table 3 (continued) 



Star 

Date 

V 

(mag) 

U-B 

(mag) 

B-V 

(nag) 

V-R 

(mag) 

V-I 

(mag) 

V-J 

(mag) 

V-H 

(mag) 

V-K 

(mag) 

M33 Var A 

Oct. 

1977 

18.55 

+0.05 

+0.84 

+1.25 

e • • 

e • e 

see 

+4.8 


Nov. 

1980 

18.45 

ess 

+1.05 

+1.10 

• • • 

+2.31 

+3.26 

+4.43 


Dec. 

1981 

18.43 

e e • 

+1.24: 

+1.14 






Aug. 

1982 

18.09 

+0.11 

+0.88 

+1.03 





M33 Var C 

Sept 

. 1976 

17.17 

-0.79 

-0.05 

see 

e • • 

see 

e • • 

e e • 


Oct. 

1977 

16.98 

-0.76 

0.00 

+0.10 

see 

e e e 

see 

see 


Nov. 

1980 

17.21 

see 

-0.01 

+0.08 

see 

e • e 

e e # 

e e e 


Dec. 

1981 

17.23 

-0.77 

+0.10 

see 

see 

e e • 

e • • 

e e • 


Aug. 

1982 

17.15 

-0.43 

+0.09 

+0.29 






Oct. 

1982 

16.48 

-0.59 

+0.09 

+0.52 






The errors in the visual photometry for most of these stars is .02 to .05 mag 


and for the infrared photometry it Is typically 0.10 mag. 


Table 4 


omSS’W 0 '* 

*OOR QUAUTY 


Summary of Luminosities and Temperatures for 
tha Hubble- Sandage Variable* 


Star 

V 

(mag) 

Date 

Ay 

(**g) 

v e 

(«•«) 

My 

_<*«£)_ 

"Ml 

<®ag) 


*v 

<■*$ 

M 

CIO' 5 K/yr' 1 ) 

M31 AF And 

16.35 

Oct. 82 

1.05 

15.30 

-8.8 

-10.6 

20000 

-9.9 

4.1 

Var A-l 

16.59 

Nov. 80 

1.05 

15.54 

-8.6 

-9.7 

15600 

-10.6 

1.4 




1.70 

14.89 

-9.2 

-11.5 

27300 

-11.3 

1.2 

AZ And 

17.79 

Dec. 81 

0.70 

17.09 

-7.0 

-8.3 

15000 

-10.2 

- 

Var 15 

17.27 

Dec. 81 

1.05 

16.22 

-7.9 

-8.4 

8300 

-8.9 

- 




1.70 

15.57 

-8.5 

-9.5 

12400 

-9.6 


M33 Var 83 

16.16 

Dec. 81 

1.10 

15.06 

-8.8 

-11.3 

29100 


3.0 


15.43 

Oct. 82 

1.10 

14.33 

-9.6 

-11.1 

30000 

-9.6 

3.0 

Ver 2 

16.36 

Nov. 80 

.97 

15.39 

-8.5 

-10.7 

22800 

-9.1 

3.5 

Var B 

16.14 

Dec. 81 

.75 

15.39 

-8.5 

-10.3 

21900 

-9.9 

- 

Var A 

18.45 

Nov. 80 

2.73 1 

15.72 

-8.2 

- 

- 

- 

- 




2 

1.02 

- 

- 

- 

- 

-9.6 

- 

Var C 

-17.2 

80/81 

0.90 

16.3 

-7.6 


_ 

-9.0 

_ 


var A - A y include* circums teller component. 
*Var A - A y interstellar only. 
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Figure 


Figure 


Figure 


Figure 


Figure 


1 The effect* of varying the extinction curve for the ultraviolet 
ere llluetrated for AF And. The energy distribution, log X f\ 
(vatts/cm ) vs. log X (microns) is shown for AF And with the 
j. 982 observations (e) corrected for ultraviolet extinction using 
the galactic curve (0) end the curve for the UfC (P). The 
extinction curves for the two galaxies only gave different 
results for wavelength shorter then about 2500&. The 2200A 
feature can be better fit with e blackbody when the IHC 
extinction curve is used. 

2 

2 The energy distribution, log X F^ (watt a/ ca. ) vs. log X (microns), 
for AF And with the 1982 observations (e). When corrected for 
the adopted extinction (0) using the galactic curve, s 20000 K 
blackbody gives a good fit to the date with a small excess in 

the infrared. 

2 

3a The energy distribution, log X F^ (watts/cm ) vs. log X (microns), 
for Var A-l with the 1980 (s) and 1981 (A) observations. A 15600 K 
blackbody is fit to the data corrected for the adopted extinction, 
open symbols. 

3b Same as Fig. 3a for Var A-l with the maximum extinction applied 

to the observations. A 27300 K blackbody la shown fit to the data 
with a small excess In the Infrared. 

2 

4a The energy distribution, log X F^ (watts/cm ) vs. log X (microns), 

for Var 83 with the Nov. 1981 observations. When corrected for 
the adopted extinction a 29100 K blackbody Is fit to the data. 
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Figure 

Figure 

Figure 

Figure 


Figure 

Figure 

Figure 


4b Saae as Fig. 4a for Var £3 vitli Cha Nov, 1982 obaervationa. 

Comparison with Fig. 4a llluetrataa the ahift in the peak flux 
for longer wavelengths as Var 83 brightened. 

2 

5 The energy distribution, log X F^ (watts/ca ) va. log X (alcrons) 
with the 1980 (e) and 1981 (4) obaervationa. When corrected for 
fcha adopted extinction (open syabols) a 22800 K blackbody gives 

a good fit to the 1980 photometry and the ultraviolet data. 

2 

6 The energy distribution, log X F ; (watts /cm ) va. log X (aicrions) 
for the 1981 observations of Var B. A 21900 K blackbody is shown 
fit to the Halted visual and ultraviolet photoaetry corrected for 
the adopted extinction. 

2 

7 The energy distribution, log X F (watts/cm ) vs. lo? X (microns) . 
for A£ And. The 1980 (e), 1981(a) and 1982 (a) visual photo 'Cwry 
is all shown. There are no ultraviolet observations. When 
corrected for the adopted extinction (open syabols) the data 

is fit reasonably well with 15000 K blackbody. 

2 

8a fhe energy distribution, log X F (watts /cm ) vs. log X (alcrons), 

for Var 15 showing tbs 1980 (e) and 1981 (4) observations. A 
8300 K blackbody la fit through the observations corrected for 
the adopted extinction (open syabols). 

8b Saae as Fig. 8a for Var 15 with the maximum extinction. 

2 

9 The energy distribution, log X F (watts/cm ' vs. log X, for 

Var A with the 1980 observations. The data is shown when 
corrected for interstellar extinction and for possible 
rlrcums teller extinction (see text). 
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Figure 10a The combined short and long wavelength ultraviolet spactra 
(1200*3000 A) of P Cyg and r\ Car In the Galaxy and of S Dor 
In the Large Magellanic Cloud (LWR only). The flux scale Is 
indicated on the left for P Cyg and on the right for n Car and 
S Dor. 

Figure 10b The observed ultraviolet spectra of Var 83 in M33 obtained at 
two different dates and of AF And in M31. The 1982 spectrum 
of Var 83 has been shifted upward by one unit and that of 
AF And by two units. These spectra have a poor signal to noise 
in the 1950 to 2350 A Interval due to the low sensitivity In 
the IITE camera in that wavelength range. 

Figure 11 The HR diagram, vs. log T^, for the most luminous stars 

(M^ £ -9.0 mag) In six Local Group galaxies. The positions 
of the H-> variables are marked with X's and by name. 
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